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Abstract 

We study the role of a very general type of flavor symmetry in controlling the strength of i?-parity violation 
in supersymmetric models. We assume that only leptons are charged under a global symmetry whose 
breaking induces lepton number (and, hence, ii-parity) violation. The charge assignments of leptons under 
this symmetry are such that the total number of independent lepton number violating couplings is reduced 
from 39 to 6. The most severe constraints on these flavor-correlated couplings arise from neutrino masses 
and mixing as well as from the non- observation of Kl — >■ /ie. We find that such a scenario predicts an almost 
vanishing smallest neutrino mass eigenvalue, allowing the upcoming generation of neutrinoless double beta 
decay experiments to shed light on the hierarchy. 



I Introduction 

i?-parity in supersymmetry is a discrete symmetry which is defined as Rp — (—1)3^+^+25'^ where B, L and S 
are the baryon number, lepton number and spin of the particle respectively (see [H El [3l |4] ) . All standard model 
(SM) particles have Rp — 1, while all supcrparticles have Rp = — 1. The assumption of i?-parity conservation in 
supersymmetric models is quite ad hoc, as this is not supported by any deep underlying principle. Historically, 
it was imposed to keep the proton stable. However, proton decay requires a simultaneous presence of B and 
L violation. Therefore, dropping all Rp violating (RpV) couplings in one go is certainly an overkill [5]. Still, 
in conventional supersymmetric theories Rp conservation is imposed primarily for the sake of convenience, as 
otherwise the number of independent parameters in the minimal supersymmetric standard model (MSSM), 
which is already very large and difficult to handle, is augmented by a set of new RpV parameters. Moreover, 
conserved Rp implies that the lightest supersymmetric particle is stable, which leads to a plausible dark matter 
candidate, and also attributes supersymmetry with a characteristic missing energy signature in colliders. On 
the other hand, if lepton number is violated, one distinct advantage is that neutrino masses can be generated 
via a perfectly rcnormalizable interaction [6l [71 [8] without the need of introducing any right-handed neutrinos. 

The most general superpotential with explicit RpV couplings is given by [3] 

Wrpv = ^i^L,Hu + \K,kUL,E'^ + K,kL^QM + iAr,,C/f Z?f I?^, (1) 

where i, j and k are the three quark and lepton generation indices. Here, Qi and Li are SU(2)-doublet quark and 
lepton superfields, respectively; Df and Ef are SU(2)-singlet superfields for down-type quarks and charged 
leptons, respectively; and, is the Higgs superfield that generates the mass of the up-type quarks. This 
introduces 48 new couplings: 3 /x^-type, 9 Aij^-type (note the antisymmetry in the first two indices), 27 A^^j.- 
type and 9 A^^-j,-type couplings (note the antisymmetry in the last two indices). Only the A" couplings are B 
violating, the rest are all L violating. Besides, new RpV soft terms appear which introduce more unknown 
parameters. We do not explicitly write down these soft terms but will mention about the relevant ones in 
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appropriate places. Dealing with so many new parameters substantially reduces the predictivity of the model. 
At this point there are two ways to proceed. One may either take one or two RpV couplings as non- vanishing at 
a time and study their implications, or, apply some suitable flavor symmetry to relate one coupling to another 
[9l [To] . We shall take the latter approach in this paper (see also fVD which does not differentiate between lepton 
doublet and down- type Higgs superfields for a different approach following a similar spirit). We show that with 
a simple flavor hypothesis we can bring down the number of totally independent RpV (more speciflcally, L 
violating) couplings to only six. These couplings induce neutrino Majorana masses, and if the neutrino mixing 
matrix is tri-bimaximal (TBM) then the number of independent L violating couplings can be further reduced 
to four, a scenario which prefers an inverse neutrino mass hierarchy. 



II A generic flavor model 

We assume that the Yukawa structure leading to the masses and mixing of quarks and charged leptons is fixed 
by some unspecified global symmetry. This symmetry also ensures baryon number conservation. There is a 
second global symmetry (X), an abelian horizontal symmetry, which is at the centre of our attention. Only 
leptons are charged under X , such that for each generation i, 

Qx{U) = -Qx{Ef). (2) 

We assume that the Qx charges of different generations are all positive. The horizontal symmetry is explicitly 
broken by a small parameter e < 1, whose charge under X is Qx(e) = — 1- If the total charge of a given 
superpotential term is n, then the term is suppressed by e". As an example, if Q = Zjsi, then the suppression 
would be e"('°°dAr) ^ 

Now we look at the consequences of Eq. ([2|) for the 48 RpV couplings of Eq. ([T]). Since i?- number is conserved, 
all the A" couplings vanish right away. Since only leptons are charged under X, it follows that Qx{LiQjD'^) = 
Qx{LiHu) = Qx{Li), and hence = A^^-^, ^ jli = ^i/^, where the supersymmetry preserving /i parameter is 
assumed to be of the same order as the supersymmetry breaking soft masses (m). Turning our attention to the 
LiLjE^ operator, we notice that when j = k, the same argument as above leads to Xijj ~ A^^^, ~ fli. Thus 39 
a priory independent L-violating couplings basically boil down to only 6: 

(— Kjk — ^ijj) J ^^123 1 ^^132, A231 . (3) 

Thanks to the flavor symmetry, the i-violating bilinear soft parameters Bi would be aligned to the corresponding 
superpotential parameters fii as well, i.e. Bi = Bijiv? ^ fii- It should be noted that when we say that two 
couplings are related, we mean that they have a common suppression factor . Indeed, there are order-one 
uncertainties in the actual coefficients of the operators, for which reason we have used a 'near-equality' sign in 
Eq. ([3]). Now we come to the relative sizes of the L- violating couplings. The suppression would depend on the 
sum of Qx charges of the associated lepton fields as a power of e. More specifically, 

/i, ~ ~ a: ~ e«-(^'), A,,fc ^ e{Q-(iO+Qx(L,)+Qx(i=;,^)} _ (4) 

Eventually, we shall provide a specific demonstration with Q = Zx^ x [9 , which means there are all together 
6 charges for the three lepton generations. 

Many RpV couplings which are not so strongly constrained may now be related by Eq. ([4]) to the ones which 
are severely constrained by experiments. The existing bounds on the individual and product couplings can be 
found in the reviews [3]. 

III Neutrino masses and mixing 

One of the high points of i?-parity violation is that it generates neutrino masses and mixing through a perfectly 
renormalizable interaction without the need of introducing any right-handed neutrino. This has already been 
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Our couplings Related to Existing limits (Sources) Refs. 

1.5 X 10-'' (M^) [m^J E~ 
Ai23 Ai23 0.03 (ffi) [Vud] [4] 

Al32 Al32 0.03 iflR) [Rr] [4] 

A231 A231 0.05 (efi) [Rr] H 

Table 1: The list of the six independent couplings and the standard couplings they are related to by the flavor symmetry 
X. The three p,i couplings are of the same order of magnitude as 36 out of 39 a priori independent RpV couplings. 
A mass of 100 GeV is assumed for the superparticles exchanged in the processes involved. These superparticles are 
indicated within first bracket right after the bounds (the weak gaugino mass AI^ and the three scalar leptons In). 
The entries in the square brackets specify the different observables from which origin the bounds originate. Here, 
Rt = r(r~ ^-F^!^T-)/r(/i- e~VeVfj). 



studied at various levels of detail [SI [71 |H]- In this work we will follow the notation of [5]. The neutrino masses, 
in the basis in which all the sneutrino vacuum expectation values vanish, can be written as 



I] ^ ^^^^J.j I g47^2 ^Qg2 p j^3 ' 547^2 ^ fj^2 

\ - 3 f^iKkk ~^ f^jKkk , \ ^ 1 fJ-iKkk + fJ-j^ikk 

+ 2^Tr^9md, z + 2^Tr^9me, 5 

^-^ Idtt-^ m ^-^ loTT^ TO 

k k 

+ -^^ik^ki — ^1 — Mtan/3 + > —^XtikXjki — r^/itan/S, 

l,k 9 l,k 

where to^^ (™ei) denote the masses of the down- type quarks (charged leptons). A comment on the approxima- 
tions made above is in order. We have denoted the squark masses by rhq and assumed them to be somewhat 
heavier than a common mass scale to assumed for the sleptons and weak gauginos/Higgsinos. This approxima- 
tion may be crude but is good enough for our order-of- magnitude estimate of the RpV couplings. In Eq. ([5]), 
the first line accounts for the tree level and one loop contributions from bilinear couplings only, the second line 
represents the one loop contributions from both bilinear and trilinear couplings, while the last line stands for 
one loop contributions from trilinear couplings only. The possibility of large left- right squark/ slepton mixing 
which may be induced by large tan/3(= Vu/vd) has been taken into account in the purely trilinear loop dy- 
namics. The tree level /^i/ij contribution generates a rank-one mass matrix and therefore yields only one mass 
eigenvalue. Since, in our case, B^, A^, are all proportional to /^i, even after including their contributions the 
rank-one nature of the mass matrix does not change. What breaks the alignment and yields more non- vanishing 
eigenvalues is the contribution from the purely trilinear loops involving Xijk{i ^ J 7^ k), since these couplings 
are not aligned with /i^. 

This leaves us with the remaining three couplings, namely, A123, A132 and A231, no two indices of which are the 
same, for generating the second mandatory and the third optional nonvanishing neutrino masses and the three 
mixing angles (two large and one small). Note that the existing bounds on Xijk with i ^ j =^ k are relatively 
less stringent ~ see Table [TJ 

Different low energy processes, especially some lepton flavor violating decays, yield important constraints on 
trilinear product couplings [12l [131 113 • the smallness of most of the couplings as shown in the first row of 

Table [U these constraints are in almost all cases easily satisfied. The bounds emerging from the nonobservation 
of /xe/e/I [131 [Hj, namely, 

>^^JkXLn < 6.7 X 10-9 mj^^ I (100 GeV)^ , (6) 

with the combinations {ijk){lmn) : (312)(312), (312)(321), (321)(312), (321)(321), play a crucial role in neutrino 
mass/mixing model building in our scenario, as we shall see later. Due to the specific inter-connections among 
RpV couplings owing to the flavor symmetry, Eq. ([6]) leads to the following limits: 

Ai32Ajj , A23iAjj < 6.7 X 10-9 ml^^ / (100 GeV)^ . (7) 
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If we set the numerical values of the couplings near their upper limits (see Table [T|) , they turn out to be large 
enough to offset the loop suppression factors. The mass matrix entries can then be written with only six RpV 
couplings as 

1 ITlnmr 



with 



ruefj. « afiifi2 + ^-^Ai23A232 — :nr^fJ'ta.nl3 + — JA213A131 ^tan^, 
OTT-^ ovr^ 

1 niniriT 1 m-u'Tie r, /c^\ 

nier « a^i^3 + — ^Ai32A323 fita.nfi+ ^^A3i2Ai2i Z^^ ^tan/J, (8) 



1 , , m^m^ 



m^^ « a^2M2 + ^-^A23iA2i3— r^A^tan^, 

1 mT-me „ 1 . , mnrrie 

^-^A23iA3i3 , ^tan/j+ — ^A32iA2i2 ^tan/^, 

OTT-^ TO^ OTT"^ 



1 ninme 
af^sfJ'S + ^-^A3i2A32i — :ri^fJ' tan 13, 



a=^ + V^ + V#?^+V|^tan/3. (9) 

k kJ^ 



With this mass matrix we try to reproduce the neutrino oscillation data, namely, the two mass-squared differ- 
ences (Ato2]^ and Amgj^) and the three mixing angles (6*12, ^23 and ^13). For simplicity we assume that all the 
phases in the neutrino mixing matrix are zero. Since neutrino oscillation analysis probes only the mass-squared 
differences and not their absolute values, we need to assume the hierarchy (normal/inverted) of the masses and 
the size of the smallest eigenvalue to fix the other two masses. There is no lower limit on the smallest neutrino 
mass eigenvalue, it can still be zero. 

We take the best fit values of the neutrino mass-squared differences from [TS]: Am^i ~ 7.59 x 10^^ eV^, 
Amii(IH) = -2.40 x lO'^ eV^ Arnii(NH) = 2.51 x lO'^ eV^ Above, NH stands for 'normal hierarchy' 
and IH stands for 'inverted hierarchy' of neutrino masses. The two mixing angles dyi and ^23 are set to their 
TBM values (using best fit values instead does not lead to significant changes) . Very recently, two long-baseline 
accelerator experiments T2K and MINOS, both probing v^^ — v^. appearance, have reported, for the first time, a 
non-zero measurement of ^13. T2K has observed 6 electron-like events against an estimated background of 1.5, 
thus discarding ^13 = at the level of 2.5cr [TB]. The MINOS experiment observes 62 electron-like events against 
an expected 49, thus disfavoring ^13 = at 1.5(7 [17]. A new global fit suggests sin^ 6*13 = 0.021(0.025) ± 0.007 
with old (new) reactor fluxes [T^ . The central value corresponds to ^13 « 9° . A second global flt can be found 
in [12]. Both indicate ^13 > with a significance of about Scr. We shall see below that in our scenario whether 
013 is vanishing or non-vanishing plays an important role in predicting whether neutrino mass hierarchy is 
normal (Am^ = m§ — 0.5(m2 -t- m^) > 0) or inverted (Arn^ < 0). Since none of the two experiments has so far 
conclusively established a nonzero value of ^13, we take both the paradigms, namely, ^13 — and (^13 ^ 0, and 
study what do they imply on the choices of RpV parameters and whether we can predict the nature of mass 
hierarchy. 



IV Tri-bimaximal mixing 

The TBM structure immediately implies that |mep| = |me^| and Iw^^l = |'71t-t-|, regardless of whether the 
lightest mass eigenvalue is vanishing or not, and also irrespective of whether the neutrino mass hierarchy is 
normal or inverted. For our couplings this can be comfortably realized by setting IA123I = IA132I and \[i2\ = ~ 
see Eq. ([5]) . This means that we can parametrize the mass matrix with four independent RpV parameters instead 
of six, which of course improves the predictivity of the model. Dropping the terms in the loop contribution 



4 



proportional to the electron mass, we obtain 

1 m,-m„ 

« TOer « a^lM2 - ■S-^Ai23M3 tan fi . (10) 

Clearly, under this situation, the absolute values for the tree-level contributions to m^;^; ^ a/^2M2j i^ht ^ 0/^2/^3 
and ~ a^^fi^ are the same. Setting all CP-violating phases to zero, the TBM mixing matrix takes the form 

m 



f/TBM= ^ ^ . (11) 




To fix the numerical values of the mass matrix from m = ?7tbm ^ diag(mi, m2, TO3) x C/tbm, all we need to 
decide is the mass hierarchy (normal or inverted) and the smallest mass eigenvalue. 

Inverted hierarchy: We first consider the case of inverted hierarchy with ma = 0. This choice additionally 
demands m^,- — — m^^. One obtains 

4.92 X 10-2 2.56 x lO"'^ -2.56 x lO""*^ 

2.56 X lO"'' 2.47 x lO'^ -2.47 x 10"^ | eV (IH, TBM, ma = 0). (12) 
-2.56 X 10-4 -2.47 x IQ-^ 2.47 x 10"^ 

By setting /i2 — — /X3 and keeping A231 ^ A123, we obtain a rough analytical solution using Eqs. (|4]) and ([8]): 



I I I I r~, , / An'^meem'^ m^^ -I- Ai2Ai23TOrTO/i tan/3/ (STr^m^) 

Ui2\ = \f^3\~ ya 'rrif^f,, A123 «W ^, fJ-i ~ r . (13) 

V W m^rrir ^jL tan p (a/^2) 

Putting fh — fi ^ 100 GeV and fhq — 300 GeV in Eq. (|13p we obtain a solution (with ^2 = — ^3, and A132 = 
— A123) 

jli = 1.9 X 10"^ , = -4.7 x 10"^ , A231 - 10"'', A123 = -3.2 x 10"* for tan/3 = 10 . (14) 

To illustrate how this coupling pattern can arise from a flavor symmetry we are providing an exemplary flavor 
group for this case. However, since this choice is not necessarily unique and our conclusions do not depend 
on the specific flavor group, we omit this exercise for the other scenarios. In this case, the required relative 
suppression can be reproduced by a Qx = Z^x Zg, family symmetry with a breaking parameter s. The necessary 
charge assignments are given by, 

QsiLi) = (2, 5) , Qb{L2) = (0, 5) , QeiLs) = (3, 2) , (15) 

which imply 

Qb{LiL2E^) = (3,0), QBiL^LsE^) = (1,2), Qb{L2L^E^) = (1,2). (16) 

These assignments lead exactly to the required suppression of the couplings with e « 0.1 as 

M2(= A^s) ^ , A^i ^ e"^ , Ai23(= A132) A231 • (17) 

In the above example, the near equality of the magnitude of the entries in the fJ- — t block is ensured by 
saturating them with the tree level contributions, while keeping the loop contributions suppressed. If, within 
the TBM framework, we now consider m3 to be slightly above zero, then m^^ ~ m-rr > K^i/jrl- To obtain 
mg = 0.001 eV with rh ^ fi ^ 100 GeV, we need fii = 1.9 x 10"^, fl2 = -4.6 x 10"^ ^3 = 4.7 x 10"^, 
A123 = —3.1 X 10"'', A132 = —3.3 X 10"*, A231 = 2.7 x lO"'^. We should note two important things: (i) These 
choices imply A231A3 = 1.3 x 10~^, which mildly overshoots the Kl — /ie bound as shown in Eq. If we 
increase further, the disagreement with the bounds deepens, (m) The 'four parameter' scenario with 
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1/^21 = ImsI Slid 1-^123 = ^132! is not compatible with a non- vanishing absolute neutrino mass scale, i.e. we cannot 
fit the data assuming these 'equalities' with ma > 0, because of the hierarchical nature of the charged lepton 
masses which appear in Eq. 

Normal hierarchy: We now consider normal hierarchy of neutrino masses. In this case the smallest mass 
eigenvalue is mi. Within the TBM structure if we keep toi = 0, it follows that m^^ = > I^V^I- The 
numerical values of the mass matrix entries are 

/ 2.90 X 10-3 2.90 X IQ-'^ -2.90 x lO^^X 
TO = 2.90 X 10-3 2.80 X 10-2 2.21 x 10"^ eV (NH, TBM, toi = 0). (18) 
\-2.90 X 10-3 2.21 X 10-2 2.80 x 10-^ / 

The couplings needed to fit these entries are /ii = —5.2 x lO-"", fl2 = 3.9 x 10-^, /is = 5.0 x 10~^, A123 — 
-4.4 X 10~3, A132 = -1.2 x 10~^, A231 = 1.0 X 10~3. Although we arc within the Kl -J> fie bound, the 
requirement rrief^ — —rripT is realized quite differently. The relative signs of the tree-level couplings invariably 
imply to*™'' k. -fTOg^°°. This difference between the experimental requirement and the tree-level contribution 
cannot be resolved, even keeping in mind that signs of the RpV couplings can be chosen at will and also 
each neutrino field can be redefined to absorb a sign. Therefore, a sign adjustment for one of the entries {e[i) 
via a large loop contribution is needed, while the loop contribution to the other one (er) becomes negligible. 
This is reflected in the large hierarchy between A123 and A132. We recall that such a sign adjustment was not 
required in the case of inverted hierarchy (TBM, TO3 = 0). If we now increase the value of toi (from zero) 
and try to fit normal hierarchy within the TBM framework, the Ki, [le bound haunts us like in the case of 
inverted hierarchy with TO3 > 0. Therefore, our most robust prediction is the tight constraint for the smallest 
mass eigenvalue. Thus, inverted hierarchy can be fit with four parameters, while normal hierarchy requires six 
parameters and a sign altering large loop correction. 



V Non-zero ^13 



In view of the recent T2K data which measures non- vanishing 6113, we study how flexible we are to accommodate 
normal or inverted hierarchies while keeping ^13 close to its central value of 9.0°. Unlike in the case of TBM 
which guarantees |TOe^| = IvTierl and to^^ — mrT, it is not possible to fit the data with 4 parameters when 
013 ^ 0. 

Inverted hierarchy: First we consider the case m3 = 0. The numerical entries of the mass matrix are given by 




4.80 X 10-2 -5.13 X 10-3 -5.63 x 10-3\ 

-5.13 X 10-3 2.53 X 10-2 -2.41 x IO-2 eV (IH, 6(13 = 9.0°, TO3 = 0) . (19) 
-5.63 X 10-3 -2.41 X 10-2 2.54 x IO-2 / 



This can be fit with jli = -1.1 x 10-^ jl2 = -4.5 x 10-^, jl^ = 4.8 x 10-^, A123 = 9.3 x 10-3, A132 = 1-1 x 10-^ 
A231 = —1.1 X 10-^. Two things are worth noting: (i) The magnitudes of A123 and A132 are separated by nearly 
three orders, while they assumed identical numerical values in the case of TBM. (ii) The tree-level contribution 
to TOe^i has the wrong sign like in the case of NH with 0i3=O. Again a large sign adjusting loop contribution 
is needed to be in agreement with the experimental data. If we now increase the value of TO3, the required 
magnitude for A231 becomes larger, and eventually beyond TO3 = 0.01 eV the Kl /ie bound overshoots. 

Normal hierarchy: For toi = 0, the mass matrix entries are given by 

/4.O6 X 10-3 8.02 X 10-3 2.29 x 10-3\ 
TO = 8.02 X 10-3 2.67 X IO-2 2.16 x IO-2 eV (NH, 6*13 = 9.0°, toi = 0) . (20) 
\2.29 X 10-3 2.16 X 10-2 2.80 x 10-2/ 

This can be reproduced with fli ^ 4.1 x lO-'^, /22 = 3.8 x 10-^, fl^ = 5.0 x 10-^, A123 = -5.3 x 10-3, 
A132 = —1.5 X 10"^, A231 = 8.3 X lO-"'. Note that A231 is small enough to satisfy the Kl — > fie bound. 
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Contrary to the case of inverted hierarchy, now no large sign-flipping correction for m^^ is needed. However, 
the difference between the values of and still leads to a hierarchy in the A-couplings. Just like in the 
previous cases, the Kl — > /^e bound begins to be relevant as soon as rrii increases to around 0.005 eV (which 
requires A231 = 1.5 x 10"*^). The main conclusion for non-zero ^13 is again that the smallest mass eigenvalue 
is required to be almost vanishing in both hierarchies. But contrary to the TBM case, now IH requires a sign 
adjustment, while NH does not. 



VI Collider signatures 

The LHC signatures of the Xijk couplings have recently been explored in [51]. In our scenario, only three 
couplings Xijk{i j 1^ k) are relatively large (10~^ — 10"''), the rest are of order 10~^. The large couplings 
are small enough to make sure that the RpV vertex is numerically relevant only at the end of a supersymmetry 
cascade when the lightest neutralino decays via a Xijk interaction, Xi ~^ l^Pv. The Xijk couplings thus give 
rise to klk or Ijlk final states plus missing energy. Depending on the numerical values of of the corresponding 
Xijk couplings the branching ratios into the lilk or Ijlk channel will scale as lAyfep. Thus both invariant mass 
distributions and number counting of the final state leptons should be a part of the search method. However, 
other decay channels like xf W^P and Xi ~^ '^^ are available due to the presence of the bilinear couplings. 
Their role has been investigated in detail in |22| . Therefore, a detailed study of neutralino decays is important 
to test this and other RpV models and differentiate between them. The non-observation of an excess in four 
lepton events at CMS and ATLAS so far indicates a somewhat heavier squark mass scale than the one we 
choose. However, scaling the slepton masses accordingly, this will not lead to any significant changes related to 
our work. 



VII Conclusions 

In this paper we have studied a generic and simple flavor model which reduces the number of independent 
couplings from 39 to 6, i.e. \jLi(i = 1,2,3), A123, A132 and A231. This results in an extremely predictive 
framework, which can reproduce the correct neutrino masses and mixings while satisfying all other low energy 
bounds. 

In its simplest realization the scenario leads to a four parameter model with exact tri-bimaximal mixing and 
prefers inverse hierarchy, for which a specific flavor model, viz Z4 x Zg, has been proposed. A non- vanishing 
mixing angle ^13 can be accommodated in a six parameter realization. 

A general prediction of all possible realizations is an almost vanishing absolute mass scale for neutrinos, i.e. an 
essentially massless lightest neutrino. This feature is tightly related to the non-observation of A'l — > /ie which 
affects many important coupling products in this framework. As a consequence, any positive signal in one of the 
upcoming neutrinoless double beta decay experiments would imply an inverted neutrino mass hierarchy, since 
for the combination of normal hierarchy and an almost vanishing absolute mass scale, the resulting |mee| is 
beyond their sensitivity. In other words, if a conclusive evidence of nonzero ^13 is established, then our scenario 
would be able to accommodate a positive signal of neutrinoless double beta decay only at the expense of large 
sign-flipping correction to one of the off-diagonal elements of the mass matrix. Moreover, the flavor structure 
proposed here can lead to specific decays of a neutralino LSP at the LHC. 



Acknowledgments: 

We thank W.Porod for valuable suggestions. This work was supported by DAAD-DST PPP Grant 
No. D/08/04933 and DST-DAAD project No. INT/DAAD/P-181/2008. G.B. acknowledges the Gambrinus 



7 



grant and hospitality at T.U. Dortmund during a part of the cohaboration. D.P. and H.P. were supported by 
DFG Grant No. PA 803/5-1, the Physics at the Terascale Helmholtz Alhancc Working Group: Neutrino masses 
and Lepton Flavor Violation at the LHC, and they acknowledge hospitality at the Saha Institute of Nuclear 
Physics, Kolkata, during a part of this collaboration. 



References 

[1] G. R. Farrar and P. Fayet, Phys. Lett. B 76 (1978) 575; S. Weinberg, Phys. Rev. D 26 (1982) 287; N. Sakai 
and T. Yanagida, Nuci. Phys. B 197 (1982) 533; C. S. Aulakh and R. N. Mohapatra, Phys. Lett. B 119 
(1982) 136. 

[2] A. Santamaria and J. W. F. Valle, Phys. Rev. Lett. 60 (1988) 397; A. Santamaria and J. W. F. Valle, Phys. 
Lett. B 195 (1987) 423; A. Santamaria and J. W. F. Valle, Phys. Rev. D 39 (1989) 1780; J. C. Romao and 
J. W. F. Valle, Nucl. Phys. B 381 (1992) 87; J. C. Romao, C. A. Santos and J. W. F. Valle, Phys. Lett. B 
288 (1992) 311. 



[3] For reviews see, G. Bhattacharyya, Nucl. Phys. Proc. Suppl. 52A (1997) 83 arXiv:hep-ph/9608415 ; 
G. Bhattacharyya, arXiv:he p-ph/9709395; H. K. Drein er, arXiv:hep-ph/9707435; M. Chemtob, Prog. 
Part. Nucl. Phys. 54 (2005) Tl~ |arXiv:hep-ph/0406029l ; R. Barbier et al, Phys. Rcpt. 420 (2005) 1 
[arXiv:hep-ph/0406039j . 

[4] For a very recent update on sing le RpV couplings, see Y. Kao and T. Takeuchi, la"rXiv:0910.4980l [hep-ph]. 

[5] Sakai and Yanagida in [T]; A. Y. Smirnov and F. Vissani, Phys. Lett. B 380 (1996) 317 
^ arXiv:hep-ph/9601387j G. Bhattacharyya and P. B. Pal, Phys. Rev. D 59 (1999) 097701 

G. Bhattacharyya and P. B. Pal, Phys. Lett. B 439 (1998) 81 



arXiv:hep-ph/9809493 



arXiv:hep-ph/9806214 



[6] L. J. HaU and M. Suzuki, Nucl. Phys. B 231 (1984) 419; M. Nowakowski and A. Pilaftsis, Nucl. 
Phys. B 461 (1996) 19 ar Xiv:hep-ph/950827l ; S. R akshit, G. Bhattacharyya and A. Raychaudhuri, 
Phys. Rev. D 59 (1999) 091701 |arXiv:hep-ph /9'811500| ; G. Bhattach aryya, H. V. Klapdor-Kleingrothaus 
and H. Pas, Phys. Lett. B 463, 77 (1999) |arXiv:hep-ph/9907432' ; D. E. Kaplan and A. E. Nelson, 
JHEP 0001 (2000) 033 arXiv:hep-ph/9901254 ; M. Hirsch, M. A. Diaz, W. P orod, J. C. Romao and 
J. W. F. Valle, Phys. Rev. D 62 (2000) 113008 [Erratum-ibid. D 65 (2002) 119901] | arXiv:hep-ph/00041 15] ; 
M. A. Diaz, M. Hirsch, W. Porod, J. C. Romao and J. W. F. Valle, Phys. Rev. D 68 (2003) 
013009 [Erratum-ibid. D 71 (2005) 059904] ,arXiv:hep-ph/0302021 ; S. Davidson and M. Losada, Phys. 
Rev. D 65 (2002) 075025 |arXiv:hep-ph/0010325|; A. Dedes, S. Rimmer and J. Rosiek, JHEP 0608 
(2006) 005 'arXiv:hep-ph/06032257X'S. Joshi^jura and M. Nowakowski, Phys. Rev. D 51 (1995) 2421 
,arXiv :hep-ph/9408224 ; R. Hempfling, Nucl. Phys. B 478 (1996) 3 ,arXiv:hep-ph/9511288 . M. Hirsch , 
H. V. Klapdor-Kleingrothaus and S. G. Kovalenko, Phys. Lett. B 398 (1997) 311 |arXiv:hep-ph/9701253] . 
Y. Grossman and H. E. Haber, Phys. Rev. Lett. 78 (1997) 3438 |arXiv:hep-ph/970242l' . Y. Grossman and 
H. E. Haber, Phys. Rev. D 59 (1999) 093008 arXiv: hep-ph/9810 536 . Neutrino mass calculation induced by 
i?-parity violation is now incorporated in SOFTSUSY code, see B. C. AUanach, C. H. Kom, M. Hanussek, 
[arXiv: 1109.37351 [hep-ph]]. 



[7] A. Abada, S. Davidson and M. Losada, Phys. Rev. D 65 (2002) 075010 [arXiv:hep-ph/0111332 ; S. David- 
son and M. Losada, JHEP 0005 (2000) 02 1 arXiv:hep-ph/0005080j A. Abada, G. Bhattacharyya and 
M. Losada, Phys. Rev. D 66 (2002) 071701 [ arXiv:hep-ph/0208009| ^ 

[8] Y. Grossman and S. Rakshit, Phys. Rev. D 69 (2004) 093002 |arXiv:he p-ph/0311310| . 

[9] T. Banks, Y. Grossman, E. Nardi and Y. Mr, Phys. Rev. D 52 (1995) 5319 |arXiv:hep-p h/9505248| . 

[10] G. Bhattac haryya, Phys. Rev. D 57 (1998) 3944 |arXiv:hep-ph/9707297); Y. Kajiyama, JHEP 0704 
(2007) 007 |arXiv:hep-ph/0702056; J. R. Elhs, S. Lola and G. G. Ross, Nucl. Phys. B 526 (1998) 115 
^arXiv:h ep-ph/9803308|. 



8 



[11] H. K. Drciner, J. Soo Kim, M. Thormeier, |arXiv:0711.43T5 ' [hep-ph]]. 

[12] G. Bhattacharyya, K. B. Chatterjee and S. Nandi, Nucl. Phys. B 831 (2010) 344 |arXiv:0911.38TIl [hep-ph]]. 

[13] H. K. Dreiner, M. Kramer and B. O'Leary, Phys. Rev. D 75, 114016 (2007) |arXiv:hep-ph/0612278| . 

[14] H. K. Dreiner, G. Polesello and M. Thormeier, Phys. Rev. D 65 (2002) 115006 |arXiv:hep-ph/0112228| . 

[15] M. C. Gonzalez-Garcia, M. Mahoni and J. Salvado, JHEP 1004 (2010) 056 |arXiv: 1001. 45241 [hep-ph]]. 

[16] K. Abe et al. [T2K CoUaboration] , IHrXmlTM^Ml [hep-ex] . 

[17] L. Whitehead [MINOS Cohaboration] , "Recent resuhs from MINOS", Joint Experimental- Theoretical Sem- 
inar (24 June 2011, fermilab, USA). 

[18] G. L. Fogh, E. Lisi, A. Marrone, A. Palazzo and A. M. Rotunno. larXiv:1106.6028 ' [hep-ph]. 

[19] T. Schwetz, M. Tortola, J. W. F. Valle, jarXiv:1108.i376l [hep-ph]]. 

[20] P. F. Harrison, D. H. Perkins, W. G. Scott, Phys. Lett. B530 (2002) 167. |hep-ph/0202074| . 

[21] N. -E. Bomark, D. Choudhmy, S. Lola, P. Osland, JHEP 1107 (2011) 070. !arXivill05jQ22] [hep-ph]] . 

[22] W. Porod, M. Hirsch, J. Romao, J. W. F. Valle, Phys. Rev. D63 (2001) 115004. |hep-ph/0011248| . 



9 



